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surface  marine  vehicles.  The  mathematical  model  is  based  on  Nomoto’s  second  order 
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provide  path  control.  A  string  of  three  vehicles  is  considered  where  each  vehicle  is  using 
the  vehicle  in  the  front  as  a  reference  point.  The  coupled  motion  stability  of  the  formation 
is  analyzed  by  linearization.  It  is  shown  that  under  the  assumed  dynamics,  guidance,  and 
control  laws,  the  stability  properties  of  the  system  decoupled  into  individual  vehicles. 
This  makes  it  possible  to  obtain  exact  analytical  results  that  can  be  used  in  design. 
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parameters  on  formation  control  and  motion  stability. 
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I.  INTRODUCTION 


A.  MOTIVATION 

Unmanned  vehicles  become  more  and  more  sophisticated  every  day.  This 
increase  in  sophistication  is  accompanied  by  an  increase  in  requirements  and  missions 
they  are  called  upon  to  complete.  As  a  result,  it  is  often  required  that  they  operate  in  a 
formation  with  specific  goals  to  accomplish.  As  the  vehicles  operate  in  formations,  it  is 
possible  that  motion  stability  may  be  an  issue.  In  many  areas,  such  as  land  vehicles,  it  has 
been  found  that  a  certain  type  of  instability,  most  notably  string  instability,  may  develop 
when  vehicles  are  travelling  in  a  formation.  Therefore,  we  need  to  investigate  the  stability 
of  motion  of  unmanned  surface  vehicles  as  they  maneuver  in  the  horizontal  plane.  We 
need  to  establish  if  they  lose  their  stability  and  under  what  conditions,  so  that  we  can 
design  or  operate  vehicles  in  a  more  efficient  and  effective  manner. 

B.  BACKGROUND 

1.  Background  and  Literature  Review 

Mathematical  models  based  on  dynamic  characteristics  are  both  necessary  and 
important  in  order  to  formulate  response  and  stability  analysis  problems  for  marine 
surface  vehicles.  Different  coordinates  systems  are  adopted  in  order  to  investigate  the 
maneuverability  of  a  marine  surface  vehicle.  Typically,  two  coordinates  are  employed  for 
modeling  the  three  degrees  of  freedom  for  vehicle  motion  control,  as  described  in  Figure 
1  [2,  7,  9].  One  is  the  earth-fixed  coordinate  system  x0y0z0,  and  the  other  is  the  body- 

fixed  coordinate  system  xyz  which  moves  together  with  the  vessel. 

There  are  several  mathematical  models  in  use  for  control  system  design  and 
analysis.  They  are  generally  linear  models  or  simplified  nonlinear  models.  Nomoto’s 
model  is  a  relatively  simple  but  effective  model  for  course  keeping  control  and  autopilot 
design  [2,  7,  9,  12].  Based  on  Nomoto’s  model,  a  state  feedback  control  law  can  be 
coupled  with  a  line  of  sight  guidance  law  to  provide  path  control  as  explained  in  [1], 
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In  the  work  reported  in  [2],  an  instability  phenomenon  arising  from  mis- 
coordination  of  guidance  and  control  laws  for  marine  surface  vehicles  moving  in  a 
platoon  was  studied.  Moreover,  the  question  of  how  it  is  possible  for  ships  traveling  in 
formation  to  exhibit  the  phenomenon  of  string  instability  was  answered. 

One  of  the  necessary  tasks  for  all  studies  involving  vehicle  maneuvering  is  an 
estimate  of  the  coefficients  in  the  equations  of  motion.  In  this  thesis  we  use  an  empirical 
Maneuvering  Prediction  Program  which  offers  two  main  options.  One  is  the  Linear 
Evaluation  which  implements  the  methods  proposed  in  [3]  for  the  assessment  of  course 
stability  and  turning  ability.  The  other  is  Turning  Prediction  which  implements  the 
multiple  linear  regression  equations  presented  in  [4]  for  the  estimation  of  turning  circle 
characteristics.  In  the  Linear  Evaluation  option,  water  depth  corrections  were  added  by 
using  regression  analysis  formulas  based  on  the  data  presented  in  [5],  In  addition,  the  trim 
corrections  used  in  the  above  program  were  presented  in  [6].  These  trim  corrections  are 
used  primarily  for  velocity  derivatives. 

2.  Thesis  Overview 

The  objectives  of  this  thesis  are  as  follows: 

First,  we  want  to  formulate  the  guidance  and  control  problem  in  the  horizontal 
plane  for  a  string  of  unmanned  vehicles  and  establish  the  conditions  for  stability.  In 
addition,  we  want  to  see  how  these  conditions  relate  to  the  conditions  for  directional 
stability  of  a  single  vehicle. 

Second,  we  want  to  investigate  the  relationship  between  such  stability  conditions 
and  fundamental  geometric  parameters  of  the  vehicle.  Such  relationships,  if  they  can  be 
derived,  will  be  very  useful  in  preliminary  design  phases. 

This  thesis  is  organized  as  follows: 

Chapter  II  contains  the  problem  formulation.  We  present  an  overview  of  the 
equations  of  motion  for  marine  vehicles  in  the  horizontal  plane  with  a  set  of  assumptions 
as  they  pertain  to  this  study.  A  reduction  of  the  order  of  the  equations  of  motion  then 
follows.  This  reduction  fonns  the  basis  for  the  development  of  control  and  guidance  laws 
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for  one  and  for  a  string  of  vehicles.  The  final  set  of  equations  that  model  the  behavior  of  a 
string  of  vehicles  in  a  string  formation  is  developed  and  presented  in  this  chapter. 

Chapter  III  presents  the  stability  analysis  of  the  final  set  of  equations  developed  in 
Chapter  2.  Stability  is  based  on  a  linearized  set  of  equations.  We  fonn  the  characteristic 
equation  of  the  system,  and  we  apply  Routh’s  criterion  in  order  to  derive  the  final  set  of 
stability  conditions.  These  are  presented  in  terms  of  operational  and  control  design 
parameters. 

Chapter  IV  presents  a  set  of  parametric  analysis  results.  We  express  the 
previously  derived  stability  conditions  in  terms  of  a  number  of  physical  parameters  of  the 
vehicles.  This  way  a  designer  can  incorporate  motion  and  fonnation  control 
considerations  into  the  early  stages  of  design  or  operation  and,  therefore,  better  match 
vehicle  capabilities  to  operational  requirements. 

Chapter  V  summarizes  the  conclusions  from  this  study  and  offers 
recommendations  for  further  research. 
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II.  PROBLEM  FORMULATION 


A.  INTRODUCTION 

A  mathematical  model  for  the  problem  of  path  keeping  in  the  horizontal  plane  is 
derived  in  this  chapter.  The  objective  of  path  keeping  in  this  context  is  to  drive  the 
surface  marine  vessel  to  follow  a  commanded  path.  Starting  with  the  equations  of  motion 
for  a  marine  surface  vessel  during  maneuvering  in  a  horizontal  plane,  we  derive  the  basic 
turning  dynamics  of  the  vessel.  By  using  these  turning  dynamics  equations,  we  suggest  a 
control  law  for  the  autopilot,  which  in  coordination  with  an  appropriate  guidance  scheme 
stabilizes  the  vessel  to  a  commanded  heading  angle  y/c  capable  of  restoring  the  vessel  to 
the  commanded  path.  Finally,  we  arrive  at  the  system  of  equations  that  describes  the 
behavior  of  three  vessels  moving  in  a  fonnation,  and  according  to  the  combined  guidance 
and  control  law. 

B.  EQUATIONS  OF  MOTION 

The  rapid  advances  in  computer  technology  have  resulted  in  successful 
applications  in  ship  engineering.  Hence,  methods  of  computer  simulation  using  the 
mathematical  models  become  increasingly  important.  In  the  design  process,  computer 
simulation  and  analysis  provide  a  convenient  tool  for  predicting  ship  maneuverability. 
Development  of  the  equations  describing  the  maneuvering  motion  is  one  of  the 
requirements  for  applying  methods  of  computer  simulation  and  analysis,  and  this 
development  process  is  presented  in  the  following  sections. 

1.  Coordinate  Systems 

The  vessel  is  considered  to  be  a  rigid  body  with  only  three  degrees  of  freedom — 
surge,  sway  and  yaw — and  which  maneuvers  in  the  horizontal  plane.  The  other  three 
degrees  of  freedom — roll,  pitch  and  heave — are  neglected  and  are  not  considered  in  this 
study  [7],  see  Figure  1.  Two  coordinate  systems  are  adopted  in  order  to  investigate  the 
maneuverability  of  a  marine  surface  vessel.  One  is  the  earth-fixed  coordinate  system 
x0y0z0 ,  and  the  other  is  the  body-fixed  coordinate  system  xyz  which  moves  together  with 
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the  vessel,  see  Figure  1.  The  x0y0  plane  and  the  xy  plane  lie  on  the  free  surface,  with  the 
x0  axis  pointing  to  the  direction  of  the  original  course  of  the  vessel,  while  the  z0  axis  and 
the  z  axis  point  upwards  vertically.  The  angle  between  the  x0  axis  and  x  axis  is  defined  as 

the  yaw  angle  y/.  When  the  maneuvering  motion  starts,  the  two  coordinate  systems 
coincide  with  each  other.  After  any  amount  of  time,  the  position  of  the  vessel  is 
determined  by  the  coordinates  x0Gand  y0Gof  the  vessel  center  of  gravity  in  the  earth- 
fixed  coordinate  system,  and  the  orientation  of  the  vessel  is  detennined  by  the  yaw  angle 
¥■ 


- >Xo 

z 


\k 

Yo 


Figure  1.  The  earth-fixed  coordinate  system  and  the  body-fixed  coordinate  system. 

After  [8]. 


The  vessel  maneuvering  motion  in  the  horizontal  plane  can  be  described  by  using 
Newton’s  second  law  and  the  yaw  rate  about  the  z  axis  which  is  defined  as  r  =  y/ .  In  the 
earth-fixed  coordinate  system,  the  equations  of  motion  are  as  follows: 

mx  o  =  X  o 
my  o  =  Y  o 
h\jjr  =  N 
6 


(1) 


where  X»  and  Y »  are  the  components  of  total  force  acting  on  the  vessel  in  the  directions 
of  x0  axis  and  v0  axis,  respectively;  N  is  the  external  moment  about  the  z  axis;  m  is  the 

vessel’s  mass;  L  is  the  moment  of  inertia  of  the  ship  about  z  axis;  Xo  and  Vo  are  the 
accelerations  in  the  directions  of  x0  axis  and  y0  axis,  respectively. 


Equations  of  motion  in  the  earth-fixed  coordinate  system  can  be  expressed  with 
respect  to  a  body-fixed  coordinate  system.  We  fix  the  origin  of  the  body-fixed  coordinate 
system  lying  on  the  center  of  gravity, 

x0G  +  x  cos  y/  -  y  sin  y/  =  x0 

y0G+xsmy/  +  ycosy/ =  y0  (2) 

z  =  z0 

Replacing  the  components  of  total  force  in  the  directions  of  x  axis  and  y  axis  by  X 
and  Y,  and  the  components  of  vessel  speed  in  the  directions  of  x  axis  and  y  axis  by  uG 

and  vG  in  Equation  (2),  we  obtain 


X0  cos  y/  +  Y0  sin  y/  -  X 
-X0  sin  y/  +  Y0  cos  y/  =  Y 


(3) 


uG  cos  y/  -  vG  sin  y/  =  x0G 
uG  sin  y/  +  vG  cos  y/  -  y0G 


(4) 


By  differentiation  of  Equation  (4)  with  respect  to  time  we  get 
uG  cos  yr  -  uGy 7  sin  yr  -  vG  sin  y/  -  vGy /  cos  y/  -  x0G 
uG  sin  y/  +  uGy /  cos  y/  +  vG  cos  y/  -  vGy /  sin  y /  =  y0G 

Substituting  Equations  (1)  and  (5)  into  Equation  (3),  we  obtain  the  equations  of 
motion  in  the  body-fixed  coordinate  system,  with  the  original  of  the  system  lying  on  the 
ship  center  of  gravity,  as 


muG  -  mvGr  =  X 
mvG  +  muGr  =  Y 

Lr  =  N  (6) 


In  practice,  it  is  more  convenient  when  the  original  of  the  body-fixed  coordinate 


system  lies  at  amidships  instead  of  the  center  of  gravity,  since  the  latter  depends  on  the 
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loading  condition  of  the  vehicle.  Assuming  that  the  ship  is  symmetrical  about  its 
longitudinal  center  plane,  the  center  of  gravity  has  the  coordinates  (xG,0,zG)  in  the 
body-fixed  coordinate  system  with  the  original  lying  at  amidships.  With  this  information 
in  mind,  the  components  of  vessel  speed  at  the  center  of  gravity,  uG  and  vG  can  be 
expressed  as  follows: 

uG  =  u 

vG=v  +  xGy  (7) 

L  =  IzG  +  mxG 

Hence,  we  obtain  the  equations  of  motion  in  the  body-fixed  coordinate  system 
with  the  original  lying  at  amidships  as  [9] 

m(u  - vr - xGr 2 )-X 
m(v  +  ur  +  xGr)  =  Y 

Izi '•  +  mxG  (v  +  ur)  =  N  ^ 

where  u  and  v  are  the  surge  and  sway  velocity,  respectively;  and  xG  is  the  x-coordinate 
of  the  center  of  mass  G. 

2.  Hydrodynamic  Forces 

In  Equation  (8),  the  components  of  total  force  and  moment  acting  on  the  surface 
marine  vessel  are  designated  as  X,  Y  and  N.  These  force  and  moment  components 
involve  the  hydrodynamic  force  and  moment  due  to  different  environmental  force  and 
moment  such  as  surrounding  water  forces,  wind  forces,  wave  forces,  rudder  and  thruster 
forces,  etc.  We  will  assume  no  other  environmental  forces. 

There  are  two  different  approaches  to  presenting  the  hydrodynamic  force  and 
moment.  One  of  them  was  proposed  by  Prof.  Martin  A.  Abkowitz,  who  was  the  Director 
of  the  MIT  Ship  Model  Towing  Tank. 

Prof.  Abkowitz  used  the  Taylor  expansion  series  to  express  the  hydrodynamic 
force  and  moment.  The  hydrodynamic  force  and  moment  are  expressed  as  the  following 
form  [7] 
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X  =  X(u,v,r,u,v,  r,  8 ) 

Y  =  Y(u,v,r,u,v,r,  8)  (9) 

N  =  N(u,  v,  r,  u,  v,  r,  8) 

These  equations  were  then  expanded  in  Taylor  series  about  the  initial  steady  state  of 
forward  motion  with  constant  speed,  u0  =U,v0  =  0 ,r0  =  0 ,u0  =  0,  v0  =  0,ro  =  0,50  =  0 .  This 
results  in  [7] 


dX ,  T77  dX  dX  dX  .  dX  .  dX  .  dX  _ 

X  -Xn  + — (u-U)  + — v  + — r  + — u  + — v  + — r  + —  8 
du  dv  dr  du  dv  dr  d8 


1 

H - 

2! 


8  .  TT.  d  d  d  .  d  .  d  .  d  ~ 

—  (u  -U)-\ - v-\ - r  H - u  H - vH - r  H - 8 

du  dv  dr  du  dv  dr  88 


~\2 


X 


15,  TT.  d  d  d  .  d  .  d  .  5  ~ 

+...H —  — (u-U)-\ - vH - r-\ - u-\ - vH - rH - 8 

n\\_du  dv  dr  du  dv  dr  88 

dY  T  dY  dY  8Y  .  dY  .  dY  .  dY  . 

7  =  K  H - (u—U)-\ - - rH - mh - vH - rH - 8 


X  +  ... 


du 


dv  dr  du  dv  dr 


1 

H - 

2! 


88 

~\2 


d  .  TT.  d  d  d  .  d  .  d  .  <3„ 

—  (u— U)-\ - r-’H - rH - u-\ - vH - rH - 8 

du  dv  dr  du  dv  dr  88 


+...  +  - 


1 


n\ 

dN 


d  d  d  d.d.d.d 

—  (i/-(7)H - H’H - r  H - u  H - v-\ - r  H - o 

du  dv  dr  du  dv  dr  88 


Y  +  ... 


dN  dN  dN  .  dN  .  dN  .  dN 


N  -Na  H - (u  —  U) H - VH - rH - it  H - VH - f  H - 8 

U  V  /  r*.  ry  .  ry  .  sy  .  ry  c* 


du 


dv  dr  du  dv 


1 

H - 

2! 


d  d  d  d.d.d.d 

—  (u  -U)-\ - v-i - r  H - w  H - vH - r  H - c> 

5n  5v  dr  du  dv  dr  88 


dr  88 

~l2 

N 


1 

+...  +  — 
n\ 


d  d  d  d.d.d.d 

—  - r;H - r  H - u  H - vH - rH - 8 

du  dv  dr  du  dv  dr  88 


(10) 


N+... 


where  u0  =U,v0  =  0,  r0  =  0,  u0  =  0,  v0  =  0,  r0  =  0, 80  =  0 . 


The  equations  of  maneuvering  motion  are  then  derived  by  substituting  Equation 
(10)  into  Equation  (8). 

The  other  approach  to  define  the  hydrodynamic  force  and  moment  was  derived  by 
the  Japanese  research  group  named  Maneuvering  Mathematical  Modeling  Group 
(MMMG).  This  approach  consists  of  the  hydrodynamic  force  and  moment  acting  on  the 
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ship  hull,  propeller,  and  rudder,  as  well  as  the  interaction  between  them.  This  expression 
can  be  described  as  in  Equation  (11). 

X  =  XH+XP+XR 

Y  —  Yh  +YP  +Yr  (11) 

N  =  Nh+Np  +  Nr 

where  the  subscripts  H,  P  and  R  refer  to  the  hull,  the  propeller  and  the  rudder, 
respectively. 

By  plugging  Equation  (11)  into  Equation  (8),  we  obtain  the  equations  of  ship 
maneuvering  motion.  This  kind  of  equation  is  called  the  MMMG  model.  If  the 
components  of  the  MMMG  model  are  expressed  in  Taylor  series  expansions  and 
truncated  to  first  order,  the  two  models  are  identical. 

C.  LINEAR  EQUATIONS  OF  MOTION 

As  mentioned  before,  the  mathematical  models  generated  by  using  equations  of  a 
marine  vessel  during  maneuvering  can  be  used  for  simulation  and  analysis.  This  helps  to 
predict  ship  maneuverability.  In  order  to  analyze  ship  maneuverability,  we  can  use  a 
simplified  or  linear  set  of  equations. 

The  force  and  moment  include  hydrodynamic  derivatives  as  coefficients. 
Assuming  that,  during  ordinary  maneuvering  motions,  the  changes  in  velocities  and 
accelerations  A u  =u-U,  Av  =  v,  A r  =  r,  A u  =u,  Av  =  v,  A r  =  r  and  the  rudder  angle  8  are 
small,  the  higher  order  terms  in  the  series  in  Equation  (10)  can  be  ignored.  Hence,  we  can 
derive  the  linear  equations  of  motion  as 

m(u  - vr  - xGr 2 )  =  X0  +  Xu (u  —U)  +  Xvv  +  X  r  +  X-ii  +  Xvv  +  X. r  +  Xs8 
m(v  +  ur + xGr)  =  Y^  +  YU  (u  -U)  +  Yvv  +  Yrr  +  Yit  +  Yv  +  Yr  +  Y§8  (12) 

Izr  +  mxG  (v  +  ur)  =  N0  +  Nu  ( u  —  U)  +  Nvv  +  Nrr  +  Nitu  +  Nvv  +  NY  +  Ns5 

Assuming  ur  =  (u-U  +  U)r  =  (it  -  U)r  +  Ur  «  Ur  ,  X()  =  0.  Y(l  =  0,  A0  =  0  and 

neglecting  the  terms  vr,xGr2  in  Equation  (12)  results  in 
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(13) 


mu  =  Xu  (u  -U)  +  Xvv  +  Xr  +  X.u  +  X.v  +  X,r  +  Xs8 
m(v  +  Ur  +  xGr)  -  Yu (u  —U)  +  Yvv  +  Yrr  +  Yu u  +  Y.v  +  Y.r  +  Ys8 
Izr  +  mxG (v  +  Ur)  -  Nu (it  —U)  +  Nvv  +  Nrr  +  N-ii  +  N-v  +  N-r  +  Ns8 

For  ease  of  understanding,  we  can  take  into  account  the  port-starboard  symmetry 
of  the  vessel.  Under  this  assumption,  many  of  the  linear  hydrodynamic  derivatives,  such 
as,  Xv,Xr,X.,X.,Xg,Yu,Y.,Nu  and  Nu  vanish.  Hence,  the  linear  equations  of  motion 
can  be  simplified  as 

(m  -  Y.  )v  -  (Y.  -  mxG )r  =  Yy  +  (F  -  mu)r  +  Ys8 
(Iz  -  N- )r  -  (N-  -  mxG  )v  =  Nvv  +  (Nr-  mxGu)r  +  Ns8  ( ,  ^ 

Non-dimensionalization  of  the  equations  of  motion  in  terms  of  water  density  p, 
vessel  length  L,  and  nominal  velocity  U  can  be  expressed  in  the  form  [9] 

(m  -  7.  )v  -  (  7  -  mxG )r  =  Yv  +  (  7  -  m)r  +  Ys8 
(4  ~N-)r- (N-  - mxG )v  =  Nvv  +  (Nr  - mxG )r  +  N s8  (]5) 

D.  AUTOPILOT  CONTROL  LAW  (NOMOTO  MODEL) 

Nomoto’s  first  order  model  is  mainly  used  to  describe  the  fundamental  turning 
dynamics  of  ship  motions  or  to  design  automatic  control  devices  such  as  autopilots.  In 
this  section,  we  will  continue  working  on  dimensionless  linear  equations  of  motion. 
Using  these  equations,  we  will  finally  obtain  Nomoto’s  first  order  model. 


1.  Turning  Dynamics 


Solving  Equation  (15)  for  v  and  v,  we  get 

(m-Y')v[Ns8-(mxGU -Ni.)r-(I:-N.)r]-(mxG-Ni)[Ys8-(mU -Yr)r-(mxG-Y.)r] 

Yv(mxG-Nv)-Nv(m-Yv) 


K  [NsS  - (mxGU  -  Nr)r  -  (/_  -  N-)r]  -  Nv  [Ys8  -  (mU  -  7,)r  -  (mxG  -  Y.)r\ 

Yv(mxG-N,)-Nv(m-Y,) 


(16) 

Differentiating  the  first  equation  of  Equation  (16)  with  respect  to  time  and  setting 
the  result  equal  to  the  second  equation,  we  obtain  the  basic  turning  dynamics  of  a  marine 
surface  vessel  [9] 
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TtT2r  +  (7)  +  T2)r  +  r  =  KS  +  KT3S 

TxT2v  +  (7)  +  T2)v  +  v  =  KvS  +  K,Tj  (17) 

where 

TT  _  (Y.  - m)(N ■  mxG )(N-  - mxG ) 

1  2  Yv  (Nr  -  mxGU)  -  Nv(Yr  -  m  U) 

T +T  _  (Yv-m)(Nr  - mxGU)  +  Yv (N-  -  I.)  (Y.-mxG)Ny+(Yr  -mU)(Nv-mx0) 

1  2  Yv{Nr  -mxGU)~ Nv(Yr - mU )  Yv{Nr -mxGU)-Nv(Yr - mU ) 

T  Ys (Ni,  - mxG )  - Ns (Y.,  -  ni) 

3  YsNv-NsYv 

T  Ys{N,-Iz)-Ns{Y,-mxG) 

4  Ys(Nr  -  mxGU)  -Ns{Yr-  m U) 

K_  YgNv  -  NsYv 

Yv(Nr-mxGU)-Nv(Yr-mU ) 

_K  _Ys(N,-mxGU)-Ng(Yr-mU) 
v  Y(Nr-mx(ZJ)-Nv(Yr-mU) 

The  first  equation  of  Equation  (17)  expresses  the  relationship  between  the  ship 
turning  rate  and  the  rudder  angle. 


We  can  further  assume  thatxG  ~  0.  T  ~  0,  Y.  «  0,  Nv  ~  0.  Mv  ~  0  .  Strictly  speaking, 

such  equations  would  be  true  for  a  fore/aft  symmetric  vehicle.  It  has  been  shown, 
however,  that  a  reasonable  degree  of  accuracy  is  maintained  even  for  vehicles  that  are  not 
symmetric.  By  using  these,  we  can  obtain 


k.-Ol 

N„ 


Ti+T2~- 


T3~~ 


T  +  T  -T  «  - 

-M  T  12  1i 


I  -  N .  m  -  Y; 


K 

m  ~  T 

K 

Y-N, 

N„ 


Y 


From  the  second  equation  of  Equation  (14)  we  get 

N.  N 


(18) 


(19) 
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By  using  T  =  TX+T2-T3  and  plugging  Equation  (18)  into  Equation  (19),  the  first 

of  Equation  (17)  can  be  further  reduced  to  a  first  order  model 

Tr  +  r  =  KS  (20) 

Equation  (20)  was  firstly  derived  by  Nomoto  by  using  the  method  of  Laplace 
Transformation.  Therefore,  this  first  order  model  is  called  the  Nomoto  model.  In  this 
model,  K  and  T  are  called  the  maneuverability  indexes,  and  they  have  explicit  relations 
with  maneuvering  characteristics. 

Equation  (18)  can  be  expressed  as  in  the  form  [1] 

1  K 

r +—r  =  — 8  or  r  =  ar +  b8  (21) 

where  the  two  parameters  are 


2.  Control  Law 

A  linear  heading  feedback  control  law  based  on  Equation  (21)  has  the  form  [1] 

S0=kx{y/-y/c)  +  k2r  (22) 

where  i//  :  is  the  commanded  heading  angle. 


By  using  \j/  =  r  ,  the  system  characteristic  equation  is  obtained  from  Equation  (21) 
and  Equation  (22).  The  system  characteristic  equation  is  as  follows 

s2  -  {a  +  bk2)s  - bkx  =  0  (23) 

The  controller  gains  can  be  computed  from  the  comparison  of  the  second  order 
system  equation  as  follows: 


s2  +2£ons  + con2  =  0 

_  a  +  2£a)n 
1  b  ’  2  b 


(24) 

(25) 


The  commanded  rudder  angle  is  given  by  [1] 

,5  =  <5„,tanh(^4 
6>  „ 


(26) 
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where  Ssat  is  the  saturation  limit  and  typically  set  at  0.4  radians,  S0  is  the  slope  of  the 
function  at  zero,  and  it  is  given  by  Equation  (22). 

E.  EQUATIONS  FOR  GUIDANCE 

The  guidance  scheme  is  described  in  [1]  and  is  basically  illustrated  in  Figure  2.  In 
Figure  2,  the  vessel  located  at  (x,y)  changes  its  direction  toward  a  target  point  D  which  is 
located  ahead  of  the  vessel  at  a  distance  d  on  the  vessel’s  nominal  path.  According  to  [1], 
pure  pursuit  guidance  is  achieved  by  commanding  a  heading  angle  y/c  equal  to  the  line  of 
sight  angle  cr 

y/c=-  tan-1  —  (27) 

d 


Figure  2.  The  Line  of  Sight  Guidance  Scheme.  After  [8] 

The  guidance  law  is  based  on  the  inertial  deviation  rate  from  the  commanded 
path  [1] 
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y  =  sin  y/ 


(28) 

It  is  worthwhile  to  note  that  guidance  law  would  be  globally  asymptotically  stable 
provided  that  the  commanded  heading  angle  y/c  is  equal  to  the  vessel’s  heading  angles  . 

Moreover,  the  commanded  heading  angle  y/c  is  a  function  of  the  vessel’s  position  and  the 

distance  d.  Therefore,  the  smaller  the  value  of  distance  d,  the  faster  the  guidance  law 
response  is.  The  autopilot  has  a  limited  reaction  time  according  to  the  specified  natural 
frequency  co„  and  damping  ratio  £.  Thus,  this  parameter  d  must  be  chosen  properly  for  the 
desired  response  to  be  achieved  [1]. 

F.  COMBINED  GUIDANCE  AND  CONTROL  LAW 

In  this  section,  the  combined  guidance  and  control  law  is  derived.  The  combined 
guidance  and  control  law  describes  the  behavior  of  a  vessel  that  has  deviated  from  the 
commanded  straight-line  path  and  attempts  to  return  by  following  a  target  point  D  (see 
Figure  2).  Then,  we  expand  this  scenario  for  three  vessels  moving  in  a  string  or  series 
formation. 

1.  One  Vessel 

Substituting  the  commanded  heading  angle  yr c  defined  in  Equation  (27)  into  the 
control  law,  Equation  (22),  we  obtain 

-1  V 

S0  -kx{y/  +  tan  — )  +  k2r 

d  (29) 

For  very  small  commanded  heading  angle  y/c.  Equation  (29)  can  be  further 
reduced  to 

Sn  =  kx  (i//  +  ^-)  +  k2r 

d  (30) 

This  is  simply  the  linearized  fonn  of  (29).  Combining  Equation  (21)  with 
Equations  (26)  and  (30)  yields 

1  v 

r  =  ar  +  bSsat  tanh( - kx  (y/  +  — )  +  k2r) 

dsat  d  (31) 

The  complete  system  of  equations  is  given  by 
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(32) 

(33) 

(34) 

(35) 

(36) 


Now  we  can  expand  the  previous  concept  for  the  case  of  three  vessels  moving  in  a 
string  or  serial  formation.  In  this  case,  the  first  vessel  attempts  to  direct  its  longitudinal 
axis  toward  a  target  point  D,  whereas  the  second  vessel  points  its  longitudinal  axis  toward 
the  first  one  and  so  on  (see  Figure  3).  In  this  consideration,  di  is  the  distance  between  the 
target  point  D  and  the  first  vessel  d3  is  the  distance  between  the  first  and  the  second 
vessel;  and  d3  is  the  distance  between  the  last  two  vessels  always  measured  on  the 
commanded  path. 
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The  system  of  equations  for  the  first  vessel  is 


Vi  =  h 

i\  =  axrx  +  bx  (kt ,  (y/x  +  %  +  k2lrx ) 
a, 

Ti=^i  (37) 

For  the  second  vessel,  they  become 
^2  =  ?2 

A  =a2r2+  b2  (kX2  (if/2  +  y\  71 )  +  k22r2 ) 

a2 

T2  =  ¥2  (38) 


r 


Figure  3.  Three  vessels  in  a  formation,  which  have  deviated  from  the  commanded 

straight  line  path  (x  axis).  After  [8] 


Finally,  the  system  of  equations  for  the  third  vessel  is 
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Vi  =r, 

1 \  =  a3V3  +  bi  (k23  (Vi  +  } ’3  ,  ~V2  )  +  V))  (39) 

a3 

j3  =¥3 

or 

x  =  Ax  (40) 

where 

x  =  \i//] ,  rt ,  v, ,  i//2 ,  r2 ,  v, ,  i//3 ,  r, ,  v3  ]  and  the  system  matrix  is  given  by 

010000000 

2  -co  2 

- CO  2  -2Cffl  — 0  0  0  0  0  0 

1  "■  dx 

100000000 

000010000 

2  2 
.  °>n  7 

A=  0  0  -co,2  -2Cxo,  — 0  0  0 

d2  "2  2  "2  d2 

000100000 

000000010 

2  2 

0  0  0  0  0  -<o  2  -2^-.co  — ^ 

d,  'h  3  "3  d, 

000000100 

(41) 
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III.  STABILITY  ANALYSIS 


A.  INTRODUCTION 

In  this  chapter,  the  problem  of  stability  for  one  vessel  is  initially  studied.  Then, 
the  problem  of  stability  for  three  vessels  moving  in  a  formation  with  constant  speed  on  a 
straight-line  commanded  path  (trivial  equilibrium  solution  characterized  by  y/=r=y=0 )  is 
analyzed.  The  solution  described  here  builds  on  and  extends  previous  work  in  [2],  All 
vessels  can  deviate  from  the  commanded  path  due  to  external  disturbances  such  as  a 
wave  or  another  external  disturbance  or  change  in  mission  requirements.  This  is 
translated  to  some  non-zero  initial  conditions.  The  system  is  prone  to  instability 
phenomena  due  to  incorrect  coordination  of  guidance  and  control  laws.  Specific  criteria 
for  stability  and  conclusions  about  the  effect  of  damping  ratio  C  and  natural  frequency  ro„ 
are  drawn. 


B. 


STABILITY  OF  ONE  VESSEL 


For  the  case  of  one  vessel,  we  have  already  referred  to  the  combined  guidance  and 
control  law,  Equation  (34),  and  the  matrix  A  is 

'0  1  0 


A  = 


V  -2$®,  - 


co„ 


1 


0 


0 


(42) 


Local  stability  properties  can  then  be  established  by  the  eigenvalues  of  matrix  A. 
The  characteristic  equation  of  A  is 

2 


CO., 


A  +  2£conA  +  oon  A  H — —  —  0 

d 


(43) 


co„ 


We  can  observe  that  if  we  neglect  the  constant  term  — —  from  Equation  (42)  (i.e., 

d 


rf->oo  or  the  guidance  law  is  eliminated),  Equation  (43)  reduces  to  Equation  (24)  (i.e., 
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the  characteristic  equation  when  control  law  stabilizes  the  vessel  to  any  commanded 
heading  angle  y/c  without  the  presence  of  the  guidance  law). 


Applying  Routh’s  criterion  to  Equation  (43)  we  get 
2  1 


co.. 


2  C(o~ - —  >  0  =>  > 

d  2$ Q)n 


d critical 


or 


(44) 


For  d  > 


1 


all  eigenvalues  of  A  have  negative  real  parts,  and  the  combined 


guidance  and  control  law  provides  stability  (i.e.,  the  vessel  follows  the  commanded  path). 


On  the  other  hand,  we  already  referred  to  the  controller  gains,  as  in  Equation  (25). 
Using  Equation  (25)  we  obtain 

(45) 


co„  =  -■> 

_  knb  +  a 


2  m 


(46) 


Substituting  Equation  (45)  and  (46)  into  Equation  (44),  dcriticaI  can  be  expressed  in 
the  form 

1 


d critical 


k2b  +  a 


(47) 


Equation  (47)  proves  that  dcritical  depends  only  on  k2  and  maneuverability 


indexes,  and  it  does  not  depend  on  kx .  In  other  words,  we  see  that  the  critical  distance  for 
stability  of  the  combined  law  is  a  function  of  the  derivative  gain  of  the  vehicle  control 
law,  an  observation  which  had  escaped  previous  studies.  We  conclude  that  the  overall 
stability  of  the  system  depends  then  on  the  speed  of  response  of  the  vehicle  to  the  rate  of 
change  of  its  commanded  heading  angle.  A  vehicle  that  has  a  higher  rate  of  change,  in 
other  words  a  higher  bandwidth  or  more  responsive  vehicle,  can  tolerate  smaller  values  of 
the  distance  d  and  still  be  stable. 


C.  STABILITIY  OF  THREE  VESSELS 

The  controller  gains  for  the  first  vessel  are  from  Equation  (25) 
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-CO. 


2 


1  0  0  0  0  0  0 

-co  2 

-2£>„  — 1  0  0  0  0  0 
1  d> 

0  0  0  0  0  0  0 

0  0  0  1  0  0  0 

2  2 
con  , 

o  2  -2^>„  — 0  0 

d2  2  2  d2 
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The  study  of  the  eigenvalues  of  matrix  A  can  reveal  the  local  stability  properties 
of  the  string.  First,  we  develop  the  characteristic  equation  of  matrix  A 
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Following  some  algebra,  we  can  show  that  this  is  the  product  of  the  following 
three  equations 

2  2 
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(53) 


Equation  (53)  is  essentially  the  same  form  as  Equation  (43),  which  is  the 
characteristic  equation  for  the  case  of  one  vessel  trying  to  move  in  the  commanded  path. 


The  characteristic  equation  of  matrix  A  then  becomes 


2  2  2 

dl  dy  ■  di  (54) 

Therefore,  the  nine  eigenvalues  of  matrix  A  are  the  roots  of  Equation  (53).  In 

other  words,  the  eigenvalues  for  the  case  of  three  vessels  independently  try  to  follow  the 

commanded  path.  This  means  that  system  stability  is  established  if  and  only  if  all  the 
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vessels  have  stability  under  the  assumption  that  they  move  independently  (not  in  a  string, 


but  separately),  as  in 
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(55) 
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IV.  PARAMETRIC  STUDIES 


A.  INTRODUCTION 

In  this  chapter,  parametric  studies  of  marine  surface  vehicle  stability  conditions 
are  introduced.  A  maneuvering  prediction  program  was  used  to  evaluate  the 
hydrodynamic  coefficients  for  a  range  of  vehicle  geometric  parameters.  Using  this 
program,  the  coefficients  in  Nomoto’s  model  was  evaluated.  Then,  we  looked  at  the 
variation  of  the  critical  parameter  for  stability  in  terms  of  vehicle  geometry. 


B.  MANEUVERING  PREDICTION  PROGRAM 


The  Maneuvering  Prediction  Program  (MPP)  we  utilized  here  was  first  developed 
to  support  the  teaching  of  conceptual  ship  design  within  the  University  of  Michigan’s 
Department  of  Naval  Architecture  and  Marine  Engineering.  This  program  applies 
methods  to  assess  the  course  stability,  turn  ability,  and  controllability  of  a  surface  marine 
vehicle.  While  applying  these  methods,  the  MPP  uses  the  following  empirical  equations 
in  Equation  (56) 
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where  B,  T,  L  are  the  maximum  beam,  the  mean  draft  and  the  length  on  waterline, 
respectively.  Y. ,  Y. ,  Nv ,  Nr ,  7 ,  F ,  Nv ,  Nr  are  the  vessel  hydrodynamic  derivatives.  These 

formulas  form  the  basis  of  the  predictions  incorporated  in  the  MPP,  although  additional, 
primarily  empirical,  adjustments  based  on  [3]  and  [4],  are  made  in  the  program. 


The  vehicle  characteristics  input  to  the  MPP  are  listed  in  Table  1. 


Vehicle  Characteristics 

Length  on  Waterline  (LWL) 

100 

Meters 

Maximum  Beam  (B) 

10 

Meters 

Draft  Forward  (TF) 

7 

Meters 

Draft  Aft  (TA) 

7 

Meters 

Block  Coefficient  on  LWL  (CB) 

0.7 

- 

Table  1.  Vessel  characteristics  used  in  MPP. 


Parametric  studies  were  performed  in  three  separate  parts.  Beam,  block 
coefficient  and  mean  draft  effects  on  critical  distances  are  analyzed.  Typical  outputs  from 
the  MPP  for  each  effect  are  included  in  the  appendix. 

C.  BEAM  EFFECTS 

In  order  to  present  the  beam  effects  on  the  critical  stability  distance,  we  kept  the 
beam  of  the  model  as  a  variable,  and  the  remaining  vehicle  characteristics  are  set  to  be 
constant  as  shown  in  Table  1.  Using  the  MPP  and  the  previous  parametric  expressions, 
we  calculated  the  T  and  K  maneuverability  indexes  for  each  value  of  different  beam 
values  assuming  that  the  surface  vehicle  was  hydrodynamically  course  stable.  Results  for 
hydrodynamically  course  unstable  vehicles  were  discarded  as  they  were  believed  to  be 
unreliable  and  could  not  be  adequately  modeled  by  the  MPP’s  parametric  formulas. 
Typical  results  for  different  beam  values  are  listed  in  Table  2. 
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Maximum  Beam  on 

Time  Constant  (T) 

Rudder  Gain  Factor  (K) 

Waterline  (B)  (meters) 

7 

0.8427 

-0.5704 

8 

1.0534 

-0.6677 

9 

1.3215 

-0.7935 

10 

1.6706 

-0.96 

11 

2.1398 

-1.1875 

12 

2.798 

-1.5117 

13 

3.7808 

-2.0028 

14 

5.3941 

-2.8191 

15 

8.5058 

-4.4104 

16 

16.9337 

-8.7544 

Table  2.  Calculated  values  of  T  and  K  maneuverability  indexes  for  different  beam 

values. 


We  already  referred  to  Nomoto’s  first  order  model  in  Equation  (20)  and  expressed 
it  in  the  form  in  Equation  (21).  Now  we  can  calculate  the  new  maneuverability  indexes,  a 
and  b,  which  have  already  been  defined  in  Equation  (21).  These  a  and  b  values  for 
different  beam  values  are  listed  in  Table  3. 
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Maximum  Beam  on 

a  (-1/T) 

b  (K/T) 

Waterline  (B)  (meters) 

7 

-1.1867 

-0.6769 

8 

-0.9493 

-0.6339 

9 

-0.7567 

-0.6004 

10 

-0.5986 

-0.5746 

11 

-0.4673 

-0.5549 

12 

-0.3574 

-0.5403 

13 

-0.2645 

-0.5297 

14 

-0.1854 

-0.5226 

15 

-0.1176 

-0.5185 

16 

-0.0590 

-0.5170 

Table  3.  Calculated  values  of  new  maneuverability  indexes  for  different  beam 

values. 


Knowing  the  a  and  b  values  allows  us  to  calculate  con  and  C  by  using  Equation 

(45)  and  Equation  (46).  The  calculations  were  made  by  using  the  following  values  for  the 
parameters: 

1 .  Controller  gain  kx  \  k]  =2 

2.  Controller  gain  k2:  k2  =  - 2.5 

Calculated  <x>  and  C  values  for  different  beam  values  are  listed  in  Table  4. 
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Maximum  Beam  on 

Waterline  (B)  (meters) 

7 

1.1635 

0.2172 

8 

1.1259 

0.2821 

9 

1.0958 

0.3396 

10 

1.0720 

0.3908 

11 

1.0535 

0.4366 

12 

1.0395 

0.4778 

13 

1.0293 

0.5148 

14 

1.0224 

0.5483 

15 

1.0183 

0.5787 

16 

1.0168 

0.6065 

Table  4.  Calculated  con  and  C  values  for  different  beam  values. 


Now  either  Equation  (44)  or  Equation  (47)  can  be  used  to  determine  the  dcritical 
values.  The  results  of  dcntical  values  for  different  beam  values  are  listed  in  Table  5. 


Maximum  Beam  on 
Waterline  (B)  (meters) 

^ critical 

7 

1.9782 

8 

1.5740 

9 

1.3433 

10 

1.1933 

11 

1.0869 

12 

1.0067 

13 

0.9435 

14 

0.8919 

15 

0.8483 

16 

0.8107 

Table  5.  Calculated  dcritical  values  for  different  beam  values. 
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The  values  in  Table  5  were  used  to  plot  the  beam  effects  on  dcritical  values.  Figure 
4  presents  dcritical  values  for  different  beam  values.  Analyzing  Figure  4  brings  us  to  our 

main  conclusion  about  beam  effects.  We  conclude  that  the  critical  distance  is  decreasing 
for  the  increasing  beam;  therefore,  wider  vehicles  can  maintain  string  formation  stability 
more  easily. 


Figure  4.  Beam  effects  on  critical  distance  values. 

D.  BLOCK  COEFFICIENT  EFFECTS 

Block  coefficient  is  the  ratio  of  the  volume  of  displacement  at  any  draft  to  the 
volume  of  a  rectangular  block  with  dimensions  equal  to  the  ship  length,  the  ship  beam 
and  the  draft.  In  order  to  analyze  the  block  coefficient  effects  on  critical  distance,  a 
similar  method  used  for  beam  effect  calculations  is  applied.  This  time  the  block 
coefficient  values  are  set  as  the  independent  variable  and  the  remaining  vehicle 
characteristics  are  set  to  be  constant  as  in  Table  1.  Using  MPP,  the  T  and  K 
maneuverability  indexes  are  calculated  for  each  value  of  different  block  coefficient 
values  for  which  marine  surface  vehicle  is  hydrodynamically  course  stable.  The  results 
for  different  CB  values  are  listed  in  Table  6. 
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Block  Coefficient  ( CB ) 

Time  Constant  (T) 

Rudder  Gain  Factor  (K) 

0.52 

0.9087 

-0.6973 

0.56 

1.0403 

-0.7438 

0.6 

1.1902 

-0.7961 

0.64 

1.3619 

-0.8553 

0.68 

1.56 

-0.9226 

0.72 

1.7902 

-1.0001 

0.76 

2.0605 

-1.0902 

0.8 

2.3814 

-1.1962 

0.84 

2.7675 

-1.3227 

0.88 

3.2397 

-1.4764 

Table  6.  Calculated  values  of  T  and  K  maneuverability  indexes  for  different  block 

coefficient  values. 


As  previously  done  in  Beam  Effects  calculations,  the  new  maneuverability 
indexes,  a  and  b,  can  be  calculated.  Calculated  a  and  b  values  for  different  block 
coefficient  values  are  listed  in  Table  7. 


Block  Coefficient  ( CB ) 

a  (-1/T) 

b  (K/T) 

0.52 

-1.1005 

-0.7673 

0.56 

-0.9613 

-0.7150 

0.6 

-0.8402 

-0.6689 

0.64 

-0.7343 

-0.6280 

0.68 

-0.6410 

-0.5914 

0.72 

-0.5586 

-0.5586 

0.76 

-0.4853 

-0.5291 

0.8 

-0.4199 

-0.5023 

0.84 

-0.3613 

-0.4779 

0.88 

-0.3087 

-0.4557 

Table  7.  Calculated  values  of  new  maneuverability  indexes  for  different  block 

coefficient  values. 
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Similarly,  <x>n  and  C  can  be  calculated  by  using  Equation  (45)  and  Equation  (46). 
For  consistency,  the  same  controller  gain  values  were  used.  Calculated  con  and  C  values 
for  different  block  coefficient  values  are  listed  in  Table  8. 


Block  Coefficient  ( CB  ) 

c 

0.52 

1.2388 

0.3301 

0.56 

1.1958 

0.3454 

0.6 

1.1566 

0.3597 

0.64 

1.1207 

0.3729 

0.68 

1.0876 

0.3850 

0.72 

1.0570 

0.3964 

0.76 

1.0287 

0.4070 

0.8 

1.0023 

0.4169 

0.84 

0.9777 

0.4263 

0.88 

0.9547 

0.4350 

Table  8.  Calculated  con  and  C  values  for  different  block  coefficient  values. 
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The  results  of  dcritical  values  for  different  block  coefficient  values  are  listed  in 
Table  9. 


Block  Coefficient  (CB) 

^ critical 

0.52 

1.2226 

0.56 

1.2104 

0.6 

1.2019 

0.64 

1.1965 

0.68 

1.1940 

0.72 

1.1933 

0.76 

1.1941 

0.8 

1.1964 

0.84 

1.1997 

0.88 

1.2039 

Table  9.  Calculated  d...,  values  for  different  block  coefficient  values. 
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Figure  5  shows  the  block  coefficient  effect  on  critical  distance.  It  can  be  seen 
from  Figure  5  that  there  is  an  optimum  block  coefficient  for  formation  stability 
properties. 
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E.  MEAN  DRAFT  EFFECTS 

In  this  section,  the  same  method  is  used  to  determine  the  maneuverability  indexes 
(T,  K,  a  and  b),  on ,  £  and  dcritical  values.  The  procedure  to  calculate  these  variables  has 
already  been  explained  in  the  previous  sections. 

The  results  of  maneuverability  indexes  T  and  K  for  different  mean  draft  values 
are  listed  in  Table  10. 


Mean  Draft  (meters) 

Time  Constant  (T) 

Rudder  Gain  Factor  (K) 

1 

3.1479 

-1.7052 

2 

4.4297 

-2.3283 

3 

4.1638 

-2.2236 

4 

3.3379 

-1.8206 

5 

2.6011 

-1.4468 

6 

2.0588 

-1.1651 

7 

1.6706 

-0.96 

8 

1.3885 

-0.809 

9 

1.178 

-0.6952 

10 

1.0166 

-0.6073 

11 

0.8896 

-0.5379 

Table  10.  Calculated  values  of  T  and  K  maneuverability  indexes  for  different  mean 

draft  values. 
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Calculated  a  and  b  values  for  different  mean  draft  values  are  listed  in  Table  11. 


Mean  Draft  (meters) 

a  (-1/T) 

b  (K/T) 

1 

-0.3177 

-0.5417 

2 

-0.2257 

-0.5256 

3 

-0.2402 

-0.5340 

4 

-0.2996 

-0.5454 

5 

-0.3844 

-0.5562 

6 

-0.4857 

-0.5659 

7 

-0.5986 

-0.5746 

8 

-0.7202 

-0.5826 

9 

-0.8489 

-0.5901 

10 

-0.9837 

-0.5974 

11 

-1.1241 

-0.6046 

Table  1 1 .  Calculated  values  of  new  maneuverability  indexes  for  different  mean  draft 

values. 


Calculated  a>n  and  C  values  for  different  mean  draft  values  are  listed  in  Table  12. 


Mean  Draft  (meters) 

1 

1.0408 

0.4979 

2 

1.0253 

0.5307 

3 

1.0335 

0.5297 

4 

1.0444 

0.5093 

5 

1.0547 

0.4769 

6 

1.0639 

0.4366 

7 

1.0720 

0.3908 

8 

1.0795 

0.3411 

9 

1.0864 

0.2883 

10 

1.0930 

0.2332 

11 

1.0997 

0.1762 

Table  12.  Calculated  con  and  C  values  for  different  mean  draft  values. 
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The  results  of  d.,.,  values  for  different  mean  draft  values  are  listed  in  Table  13. 


Table  13.  Calculated  values  for  different  mean  draft  values. 


Figure  6  represents  the  block  coefficient  effects  on  critical  distance.  Focusing  on 
Figure  6,  we  can  say  that  a  shallow  draft  is  preferable  although  an  optimum  may  exist  for 
formation  stability  properties. 
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Figure  6.  Mean  draft  effects  on  critical  distance  values. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

This  thesis  addressed  the  problem  of  coordinated  motion  control  and  the  stability 
loss  of  surface  marine  vehicles.  The  problem  was  decoupled  into  the  motion  stability 
control  task  of  making  each  marine  surface  vehicle  follow  a  target  along  its  path. 

The  main  work  and  results  are  summarized  as  follows: 

1.  A  mathematical  model  was  derived  which  is  based  on  Nomoto’s  second 
order  model  and  captures  the  fundamental  dynamics  of  turning  on  the 
horizontal  plane  with  no  side  slip. 

2.  A  state  feedback  control  law  was  coupled  with  a  line  of  sight  guidance  law 
to  provide  path  control.  A  string  of  three  vehicles  was  considered  where 
each  vehicle  is  using  the  vehicle  in  the  front  as  a  reference  point. 

3.  The  coupled  motion  stability  of  the  fonnation  was  analyzed  by 
linearization.  It  was  shown  that  under  the  assumed  dynamics,  guidance, 
and  control  laws,  the  stability  properties  of  the  system  decouple  into 
individual  vehicles.  This  makes  it  possible  to  obtain  exact  analytical 
results  that  can  be  used  in  design. 

4.  Parametric  runs  and  sensitivity  analysis  studies  revealed  the  effects  of 
main  vehicle  geometric  parameters  on  formation  control  and  motion 
stability. 

5.  It  was  established  that  the  critical  stability  coefficient  is  a  function  of  the 
derivative  gain  of  the  control  law. 

6.  The  critical  stability  coefficient  depends  as  follows  on  vehicle  geometry: 

a.  It  is  decreasing  for  increasing  vehicle  beam. 

b.  It  is  in  general  increasing  for  increasing  draft. 

c.  It  has  an  optimum  in  terms  of  the  block  coefficient. 
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B.  RECOMMENDATIONS 

The  study  presented  in  this  thesis  can  be  extended  in  different  areas,  briefly 
explained  below. 

The  main  recommendations  for  future  research  are  summarized  as  follows: 

1 .  Analyze  the  nonlinear  behavior  of  the  system  in  order  to  see  if  there  are 
higher  order  effects  that  could  not  be  captured  by  linearization. 

2.  Expand  the  range  of  parametric  studies  to  additional  variables  such  as 
forward  speed  and  develop  a  set  of  design  recommendations. 

3.  Finally,  incorporate  the  longitudinal  equations  of  motion  by  allowing  the 
distances  di,  d2,  and  d3  to  be  functions  of  time. 
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University  of  Michigan 

Department  of  Naval  Architecture  ana  Marine 

Engineering 

Maneuvering  Prediction  Program  (MPP-1.3)  Dy  M.G.  Parsons 

References:  Clarke. D. .  'Sealing, P. .  andHlne.G., 

"The  Application  of  Manoeuvring  Criteria  In  Hull 

Design  using  Linear  Theory. " 
Lyster,  C.,  ana  Knlgnts,  H.  L. 

Trans. 

.  RINA.  1983 

"Prediction  Equations  for  Snips"  Turning  Circles," 

Trans.  NECIES.  1978-1979 

Bun  identification:  Beam  Effects 

input  verification: 

Length  or  tfaterllne  LML  (m) 

- 

100.00 

Maximum  Beam  on  LML  (m) 

- 

12.00 

Mean  Drart  (■) 

- 

7.00 

Draft  Forward  (m) 

- 

7.00 

Draft  Aft  (mi 

- 

7.00 

Block  Coerriclent  on  LML  CB 

- 

0.7000 

Molded  Volume  (nT3) 

- 

5880.00 

Center  or  Gravity  LOG  (%LML;  ♦  FVd) 

0.0000 

Center  of  Gravity  LCG  (m  frocs  FP) 

- 

50.00 

Midships  to  Rudder  CE  XR  ( %LML ,•  ♦  Art) 

49.0000 

Rudder  Center  of  Effort  XR  (b  from  FP) 

- 

99.00 

initial  snip  speed  (knots) 

15.00 

initial  snip  speed  (b/s> 

- 

7.7166 

Mater  Type 

SaltulSC 

Mater  Density  (kg/m* 3) 

1025.87 

Kinematic  Viscosity  (m*2/s) 

- 

0.11B831E-05 

Yaw  Radius  of  Gyration  K33/LML 

- 

0.2500 

Mater  Depth  to  Ship  Draft  Ratio  H/T 

- 

1000.00 

Steering  Gear  Tlee  Constant  <s) 

- 

2.50 

Total  Rudder  Area  -  Fraction  of  LWL"T 

- 

0.0136 

Number  of  Propellers 

- 

1 

Type  of  Single  Screw  Stern 

- 

Closed 

Sunmerged  Bow  Area  -  Fraction  of  LML«T 

0.0000 

1 
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Project  Nar.e ,  Beam  Streets 


Sat  jun  OB  13:47:40  2013 


University  or  Michigan 

Department  or  Naval  Architecture  ana  Marine  Engineering 
Maneuvering  Prealctlon  Program  (MPP-1.3)  by  M.G.  Parsons 


Linear  Maneuvering  Criteria  Option 

Rererence:  Clarice, D..  Geailng.P.,  ana  Hlne.G., 

•The  Application  or  Manoeuvring  Criteria  In  Hull 
Design  using  Linear  Theory,"  Trans.  RINA,  1983 


Run  iaentltlcatlon:  Beam  Errects 


Linear  Maneuvering  Derivatives 


Kondlaensional  Mass 

V 

prime 

. 

0 . 011760 

Nondlmenslonal  Mass  Moment 

I 

SUE 

zz 

- 

0.000735 

Sway  Velocity  Derivative 

Y 

SUE 

V 

-0 . 023640 

Sway  Acceleration  Derivative 

Y 

SUE 

V 

dot 

- 

-0.017219 

Yaw  Velocity  Derivative 

N 

SUE 

V 

m 

-0.009663 

Yaw  Acceleration  Derivative 

N 

SUE 

V 

dot 

• 

-0.000950 

Sway  velocity  Derivative 

Y 

SUE 

r 

m 

0.006164 

Sway  Acceleration  Derivative 

Y 

SUE 

r 

dot 

- 

-0 .001066 

Yaw  Velocity  Derivative 

N 

SUE 

r 

• 

-0.004049 

Yaw  Acceleration  Derivative 

N 

SUE 

r 

dot 

- 

-0.000987 

Sway  Rudder  Derivative 

Y 

SUE 

delta 

• 

0.002856 

Yaw  Rudder  Derivative 

N 

SUE 

delta 

* 

-0.001399 

Time  Constanta  and  Gains  for 

Notnoto' 

■s 

Equation 

Dominant  Ship  Time  Constant 

T1  prime 

• 

3.1175 

Ship  Time  Constant 

T2  prime 

• 

0.3669 

Numerator  Time  Constant 

T3  prime 

• 

0.7064 

Numerator  Time  Constant 

T4  prime 

- 

0.3321 

1st  order  Eqn.  Time  Constant 

T  prime 

2.7980 

Rudder  Gain  Pactor 

K  prime 

• 

-1.5117 

Rudder  Gain  Pactor 

K  sue  v  prime  - 

0.4787 

Steering  Gear  Time  Constant 

TE  prime 

- 

0.1929 

Evaluation  or  Turning  Ability  and  Stability 

Inverse  Time  Constant 

1/ IT  prime 

• 

0.3574 

Inverse  Gain  Pactor 

1/ j K  prime 

- 

0.6615 

Clarke’s  Turning  Index 

P  - 

0.3065 

Linear  Dynamic  Stability  Criterion 

C  -  0. 

.0000405 

Vessel  Is  hyOrodynamlcally  open  loop  course  staple 

Closed  Loop  Phase  Margin  with  Steering  Engine  •  28.5803  degrees 


2 
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Project  N’are :  Beam  Effects 


Sat  Jun  09  13:47:40  2013 


University  of  Michigan 

Department  of  Naval  Architecture  ana  Karine  Engineering 
Maneuvering  Prediction  Program  (MPP-1.3)  by  M.G.  Parsons 

Turning  Prediction  Option  •** 


Reference:  Lyster,  C.(  and  Knights,  H.  L. , 

■Prediction  Equations  tor  Ships*  Turning  Circle", 
Trans.  NECIES,  1978-1979 


Run  Identification:  Beam  Effects 

Approach  Speed 

- 

15.00 

Knot  8 

Rudder  Angle 

- 

10.00 

degrees 

Steady  Turning  Diameter 

- 

758.41 

meters 

Tactical  Diameter 

• 

791.10 

meters 

Advance 

• 

540.97 

meters 

Transrer 

• 

3B6.8C 

meters 

Steady  speed  in  Turn 

- 

11.19 

Knots 

3 
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Project  N’ace :  Beam  Effects 


Sat  Jun  09  13:47:40  2013 


Department 

Maneuvering 


University  ot  Michigan 
of  Naval  Architecture  ana  Karine  Engineering 

Preaictlon  Program  (MPP-1.3)  by  M.G.  Parsons 


•••State  variable  Equations  Option  ••• 


Run  Identification:  Beam  Effects 


[pal 
a/at[  r  ] 
[  v  i 


[01  01  tpsl] 

[  0  a22  a23  1 [  r  ]  + 
[  0  a32  a33  1  [  v  1 


[0  1 

[02!  delta 
(b3i 


10  0  1 

*  [garma  21  gamma  22] (y  external] 

[gamma  31  gamma  32] [N  external] 


Open  Loop 
open  Loop 
Open  Loop 
Open  Loop 


Dynamics  Matrix  Coeff. 
Dynamics  Matrix  Coeff. 
Dynamics  Matrix  Coeff. 
Dynamics  Matrix  Coeff. 


a22  - 

-2.29168 

a23  - 

S.3BB56 

a32  - 

0.17521 

333  - 

-0.61337 

Control  Distribution  Matrix  Coeff. 
Control  Distribution  Matrix  Coeff. 


b2  •  -0.99527 
b3  -  -0.13190 


EXt. 

Force 

Dlstr. 

Matrix 

Coerr. 

gaana 

21  - 

-19.43766 

Ext. 

Force 

Dlstr. 

Matrix 

coerr. 

gasma 

22  - 

592.91730 

Ext. 

Force 

Dlstr. 

Matrix 

coerr. 

gaana 

31  - 

-35.23790 

EXt. 

Force 

Dlstr. 

Matrix 

coerr. 

ga=ma 

32  - 

22.26B46 

4 


Figure  7.  One  of  the  outputs  of  the  MPP  for  beam  effects. 


44 


Project  Name:  Bloc*  Coerriclent  Effects  Sat  Jun  08  14:02:24  2013 


Onlverslty  or  Michigan 

Department  ot  Naval  Architecture  ana  Marine 

Engineering 

Maneuvering  Preaictlon  Program  (MPP-1.3)  by  M.Q.  Parsons 

References:  Clarke, D..  Geaimg.P.,  anaHlne.G., 

■The  Application  ot  Manoeuvring  Criteria  In  Hull 

Design  using  Linear  Theory. ’ 
Lyster,  C.,  ana  Knights,  H.  L. 

Trans 

.  RINA,  1983 

•Preaictlon  Equations  tor  Ships'  Turning  Circles,* 

Trans.  NECIES,  1978-1979 

Run  lOentltlcatloni  Block  Coerriclent  Effects 

Input  Verl ricat Ion . 

Length  of  Waterline  LWL  (m) 

- 

100.00 

Maximum  Beam  on  LWL  (m) 

- 

10.00 

Mean  Draft  (■) 

- 

7.00 

Drart  Forward  (m) 

- 

7.00 

Draft  Aft  (ml 

- 

7.00 

Block  coefficient  on  lwl  CB 

- 

0.8000 

Molded  volume  (m*3) 

- 

5600.00 

Center  of  Gravity  LOG  (»lwl;  ♦  Fwd) 

0.0000 

Center  of  Gravity  LOG  (m  rroca  FP) 

- 

50.00 

Midships  to  Rudder  CE  XR  ( %LKL ;  ♦  Aft) 

• 

49.0000 

Ruaaer  Center  of  Effort  KR  (■  from  FP) 

- 

99.00 

initial  Ship  Speea  (knots) 

• 

15.00 

initial  Ship  Speea  (m/s) 

7.7166 

Water  Type 

SaltulSC 

water  Density  (kg/m*3) 

- 

1025.87 

Kinematic  Viscosity  (m*2/s) 

- 

0 . 11B831E-05 

Yaw  Raaius  of  Gyration  K33/LWL 

- 

0.2500 

Water  Depth  to  Ship  Draft  Ratio  H/T 

- 

1000.00 

Steering  Gear  Time  Constant  (s) 

- 

2.50 

Total  Ruaaer  Area  -  Fraction  of  LWL*T 

- 

0.0125 

Number  of  Propellers 

- 

1 

Type  of  SlDgle  Screw  Stern 

- 

Closed 

Submerged  Bow  Area  -  Fraction  ot  LWL*T 

0.0000 

1 
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Project  Maces  BlocK  Coefficient  Effects 


Sat  jun  OB  14:02:24  2013 


University  or  Michigan 

Department  or  Naval  Architecture  ana  Marine  Engineering 
Maneuvering  Prealctlon  Program  (MPP-1.3)  by  M.G.  Parsons 


Linear  Maneuvering  Criteria  Option 

Reference:  ClarKe.D.,  Sealing. P.,  ana  Blne.G., 

•Tne  Application  or  Manoeuvring  Criteria  In  Bull 
Design  using  Linear  Theory,’  Trans.  RINA,  1983 


Run  laentiricatlon:  Block  coefficient  Effects 


Linear  Maneuvering  Derivatives 


Konaioensional  Mass 

V 

prime 

. 

0.011200 

Nonaimenslonal  Mass  Moment 

I 

SUE 

zz 

- 

0.000700 

Sway  Velocity  Derivative 

Y 

SUE 

V 

-0 . 02321B 

Sway  Acceleration  Derivative 

Y 

SUE 

V 

dot 

- 

-0.017424 

Yaw  Velocity  Derivative 

N 

SUE 

V 

m 

-0.009897 

Yaw  Acceleration  Derivative 

N 

SUE 

V 

dot 

• 

-0.000792 

Sway  velocity  Derivative 

Y 

SUE 

r 

m 

0.006455 

Sway  Acceleration  Derivative 

Y 

SUE 

r 

dot 

- 

-0 . 00092  B 

Yaw  Velocity  Derivative 

N 

SUE 

r 

• 

-0.004033 

Yaw  Acceleration  Derivative 

N 

SUE 

r 

dot 

- 

-0.001074 

Sway  Ruaoer  Derivative 

Y 

SUE 

delta 

0.002625 

Yaw  Ruaaer  Derivative 

N 

SUE 

delta 

* 

-0.0012B6 

Time  Constanta  ana  Gains  for 

Notnoto' 

■s 

Equation 

Dominant  Ship  Time  Constant 

T1  prime 

• 

2.6776 

Ship  Time  Constant 

T2  prime 

• 

0.4003 

Numerator  Time  Constant 

T3  prime 

• 

0.6965 

Numerator  Time  Constant 

T4  prime 

- 

0.3505 

1st  Oraer  Eqn.  Time  Constant 

T  prime 

2 .3814 

Ruaoer  Gain  Pactor 

K  prime 

• 

-1.1962 

Ruaaer  Gain  Pactor 

K  sue  v  prime  - 

0.3575 

Steering  Gear  Time  Constant 

TE  prime 

- 

0.1929 

Evaluation  of  Turning  Ability  ana  Stability 

Inverse  Time  Constant 

1/ IT  prime 

• 

0.4199 

Inverse  Gain  Pactor 

1/|K  prime 

- 

0.8360 

Clarke’s  Turning  InOex 

P  - 

0.2706 

Linear  Dynamic  Stability  Criterion 

C  -  0. 

.0000467 

Vessel  Is  hyOrodynamlcally  open  loop  course  staple 

Close-3  Loop  Phase  Margin  with  Steering  Engine  -  33.0114  aegrees 
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Project  Name ;  BlocK  Coerriclent  Er recta 


Sat  Jun  09  14:02:24  2013 


Unlveralty  or  Michigan 

Department  or  Naval  Architecture  ana  Karine  Engineering 
Maneuvering  Prediction  Program  (MPP-1.3)  by  M.G.  Paraons 

Turning  Prediction  Option  •** 


Reterence:  Lyster,  C.(  and  Knlghta,  H.  L. , 

■Prediction  Equations  tor  Ships’  Turning  Circle’, 
Trans.  NECIES,  1978-1979 


Run  Identlllcatlon:  Block  Coerricient  Streets 


Approach  Speed 

- 

15.00 

Knot  8 

Rudder  Angle 

- 

10.00 

degrees 

Steady  Turning  Diameter 

. 

585.81 

meters 

Tactical  Diameter 

• 

€34.03 

meters 

Advance 

• 

459. 4€ 

meters 

Transrer 

• 

308.74 

meters 

Steady  speed  in  Turn 

- 

9.41 

Knots 
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Project  Nare :  BlocK  Coefficient  Effects 


Sat  Jun  09  14:02:24  2013 


Department 

Maneuvering 


University  of  Michigan 
of  Naval  Architecture  ana  Marine  Engineering 

Preaictlon  Program  (MPP-1.3)  by  M.G.  Parsons 


•••State  variable  Equations  Option  ••• 


Run  IOentlf lcatlon-.  Block  Coefficient  Effects 


ipsi: 

a/at[  r  i 
[  v  i 


[01  01  tpsi]  [0  1 

[  0  a22  a23  1 [  r  1  •  [b2!  aelta 
[  0  a32  a33  1  [  v  1  [D3] 


(0  0  1 

•  [gasana  21  gamma  22] (Y  external] 

[gamma  31  gamma  32] [N  external] 


Open  Loop 
open  Loop 
Open  Loop 
Open  Loop 


Dynamics  Matrix  Coeff. 
Dynamics  Matrix  Coeff. 
Dynamics  Matrix  Coeff. 
Dynamics  Matrix  Coeff. 


a22  - 

-2.23190 

a23  - 

5.29389 

a32  - 

0.14828 

a33  - 

-0.63959 

Control  Distribution  Matrix  Coeff.  b2  -  -0.77727 

Control  Distribution  Matrix  Coeff.  b3  -  -0.11690 


EXt. 

Force 

Dlstr. 

Matrix 

coerr. 

gaana 

21  - 

-15.82067 

Ext. 

Force 

Dlstr. 

Matrix 

coerr. 

gasma 

22  - 

572.00311 

Ext. 

Force 

Dlstr. 

Matrix 

coerr. 

gaana 

31  - 

-35.44898 

EXt. 

Force 

Dlstr. 

Matrix 

coerr. 

gacma 

32  - 

18.53906 

4 


Figure  8.  One  of  the  outputs  of  the  MPP  for  block  coefficient  effects. 
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Project  Blase:  Mean  Draft  Effects  Sat  Jun  OB  14:09:33  2013 


Onlverslty  or  Michigan 

Department  ot  Naval  Architecture  ana  Marine 

Engineering 

Maneuvering  Preaictlon  Program  (MPP-1.3)  by  M.Q.  Parsons 

References:  Clarke, D..  Geaimg.P.,  anaHlne.G., 

■The  Application  ot  Manoeuvring  Criteria  In  Hull 

Design  using  Linear  Theory. ■ 
Lyster,  C. ,  ana  Knights,  H.  L. 

Trans 

.  RINA,  1983 

•Preaictlon  Equations  tor  Ships'  Turning  Circles,* 

Trans.  NECIES,  1978-1979 

Run  lOentltlcatloni  Mean  Draft  Errects 

Input  Verl ricat Ion . 

Length  of  Waterline  DHL  (m) 

- 

100.00 

Maximum  Beam  on  LWL  (m) 

- 

10.00 

Mean  Draft  (■) 

- 

5.00 

Drart  Forward  (m) 

- 

5.00 

Draft  Aft  (ml 

- 

5.00 

Block  coefficient  on  lwl  CB 

- 

0.7000 

Molded  volume  (m*3) 

- 

3500.00 

Center  of  Gravity  LOG  (*lwl;  ♦  Fwd) 

0.0000 

Center  of  Gravity  LOG  (m  rroca  FP) 

- 

50.00 

Midships  to  Rudder  CE  XR  ( %LKL ;  ♦  Aft) 

• 

49.0000 

Ruaaer  Center  of  Effort  KR  (■  from  FP) 

- 

99.00 

initial  Ship  Speea  (knots) 

• 

15.00 

initial  Ship  Speea  (m/s) 

7.7166 

Water  Type 

S altulSC 

water  Density  (kg/m*3) 

- 

1025.87 

Kinematic  Viscosity  (m*2/s) 

- 

0 . 11B831E-05 

Yaw  Raaius  ot  Gyration  K33/LWL 

- 

0.2500 

Water  Depth  to  Ship  Draft  Ratio  H/T 

- 

1000.00 

Steering  Gear  Time  Constant  (s) 

- 

2.50 

Total  Ruaaer  Area  -  Fraction  of  LWL*T 

- 

0.0125 

Number  of  Propellers 

- 

1 

Type  of  SlDgle  Screw  Stern 

- 

Closed 

Submerged  Bow  Area  -  Fraction  ot  LWL*T 

0.0000 

1 
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Project  Mare ,  Mean  Draft  Effects 


Sat  Jun  OB  14:09:33  2013 


University  of  Michigan 

Department  or  Naval  Architecture  ana  Marine  Engineering 
Maneuvering  Prealctlon  Program  (MPP-1.3)  by  M.G.  Parsons 


Linear  Maneuvering  Criteria  Option 

Reference:  ClarKe.D.,  (Sealing, P.,  ana  Bine, G., 

•The  Application  or  Manoeuvring  Criteria  In  Bull 
Design  using  Linear  Theory,’  Trans.  RINA,  1983 


Run  laentlflcatlon:  Mean  Draft  Effects 


Linear  Maneuvering  Derivatives 


Kondlaecsional  Mass 

V 

prime 

. 

0.007000 

Nondlmenslonal  Mass  Moment 

I 

SUE 

zz 

- 

0.000437 

Sway  Velocity  Derivative 

Y 

SUE 

V 

-0 . 012615 

Sway  Acceleration  Derivative 

Y 

SUE 

V 

dot 

- 

-0.009213 

Yaw  Velocity  Derivative 

N 

SUE 

V 

-0.004594 

Yaw  Acceleration  Derivative 

N 

SUE 

V 

dot 

- 

-0.000220 

Sway  velocity  Derivative 

Y 

SUE 

r 

• 

0.003731 

Sway  Acceleration  Derivative 

Y 

SUE 

r 

dot 

- 

-0.000423 

Yaw  Velocity  Derivative 

N 

SUE 

r 

-0.002271 

Yaw  Acceleration  Derivative 

N 

SUE 

r 

dot 

- 

-0.0005B2 

Sway  Rudder  Derivative 

Y 

SUE 

delta 

0 . 001875 

Yaw  Rudder  Derivative 

N 

SUE 

delta 

* 

-0.000919 

Time  Constanta  and  Gains  for 

Notnoto' 

■s 

Equation 

Dominant  Ship  Time  Constant 

T1  prime 

• 

2.9575 

Ship  Time  Constant 

T2  prime 

• 

0.3945 

Numerator  Time  Constant 

T3  prime 

• 

0.7509 

Numerator  Time  Constant 

T4  prime 

- 

0. 316B 

1st  order  Eqn.  Time  Constant 

T  prime 

2.6011 

Rudder  Gain  Pactor 

K  prime 

• 

-1.446B 

Rudder  Gain  Pactor 

K  sue  v  prime  - 

0.5154 

Steering  Gear  Time  Constant 

TE  prime 

- 

0.1929 

Evaluation  of  Turning  Ability  and  Stability 

Inverse  Time  Constant 

1/ IT  prime 

0.3B45 

Inverse  Gain  Factor 

1/ j K  prime 

- 

0.6912 

Clarke’s  Turning  Index 

P  - 

0.316B 

Linear  Dynamic  Stability  Criterion 

C  -  0. 

.0000141 

Vessel  Is  hyOrodynamlcally  open  loop  course  stable 

Closed  Loop  Phase  Margin  with  Steering  Engine  -  30.7841  degrees 
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Project  Nane :  Mean  Draft  Effects 


Sat  Jun  09  14:09:33  2013 


University  of  Michigan 

Department  of  Naval  Architecture  ana  Karine  Engineering 
Maneuvering  Prediction  Program  (MPP-1.3)  by  M.G.  Parsons 

Turning  Prediction  Option  •** 


Reference:  Lyster,  C.(  and  Knights,  H.  L. , 

■Prediction  Equations  tor  Ships*  Turning  Circle*, 
Trans.  NECIES,  1978-1979 


Run  Identification:  Kean  Draft  Effects 


Approach  Speed 

- 

15.00 

Knot  8 

Rudder  Angle 

- 

10.00 

degrees 

Steady  Turning  Diameter 

. 

784.83 

meters 

Tactical  Diameter 

• 

815.14 

meters 

Advance 

• 

553.45 

meters 

Transrer 

• 

398.75 

meters 

Steady  speed  in  Turn 

11.46 

Knots 
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Project  Naze:  Mean  Draft  Effects 


Sat  Jun  09  14:09:33  2013 


Department 

Maneuvering 


University  of  Mlcnigan 
of  Naval  Architecture  ana  Marine  Engineering 

Preaictlon  Program  (MPP-1.3)  by  M.G.  Parsons 


•••State  variable  Equations  Option  ••• 


Run  Identification:  Mean  Draft  Effects 
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•  [b2]  delta 
Ib3i 
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open  Loop 
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Figure  9.  One  of  the  outputs  of  MPP  for  mean  draft  effects. 


52 


LIST  OF  REFERENCES 


[1]  F.  A.  Papoulias,”  On  the  nonlinear  dynamics  of  pursuit  guidance  for  marine 
vehicles,”/,  of  Ship  Research,  vol.37,  no. 4,  pp. 342-353,  1993. 

[2]  C.  K.  Angelopoulos,  “String  stability  of  multiple  surface  marine  vessels,”  M.S. 
thesis,  Mechanical  and  Aerospace  Engineering  Department,  NPS,  Monterey, 
California,  2011. 

[3]  D.  Clarke,  P.  Gedling,  and  G.  Hine,  “The  application  of  manuevering  criteria  in 
hull  design  using  linear  theory,”  Transactions  RINA,  vol.  124,  1982. 

[4]  C.  A.  Lyster,  and  H.  L.  Knights,  “Prediction  equations  for  ships'  turning  circles,” 
Transactions  NECIES,  1978-1979. 

[5]  M.  Fugino,  “Maneuverability  in  restricted  waters:  State  of  the  art,”  The  University 
of  Michigan,  Report  No.  184,  Aug.,  1976. 

[6]  S.  Inoue,  M.  Hirano,  and  K.  Kijima,  “Hydrodynamic  derivatives  on  ship 
manuevering,”  International  Shipbuilding  Progress,  vol.  28,  no.  321,  1981. 

[7]  Z.  Zaoijan,  “Ship  maneuvering  and  seakeeping,”  Lecture  Notes,  Shanghai  Jiao 
Tong  University,  2006. 

[8]  C.  Daskam,  “DDG-2025:  Combating  the  FIAC  threat,”  Project  Report,  NPS, 
Monterey  California,  2012. 

[9]  F.  A.  Papoulias,  “Dynamics  of  marine  vehicles,”  Lecture  Notes,  NPS,  Monterey 
California,  1993. 

[10]  F.  A.  Papoulias,  “Cross  track  error  and  proportional  turning  rate  guidance  of 
marine  vehicles,”  J.  of  Ship  Research,  vol. 38,  p.  123,  1994. 

[11]  F.  A.  Papoulias,  “Bifurcation  analysis  of  line  of  sight  vehicle  guidance  using 
sliding  modes,”  Int.  J.  of  Bifurcation  and  Chaos,  vol.  1,  p.4,  1991. 

[12]  F.  A.  Papoulias,  “Stability  considerations  of  guidance  and  control  laws  for 
autonomous  underwater  vehicles  in  the  horizontal  plane,”  in  Proceedings  7th 
International  Symposium  on  Unmanned  Underwater  Submersible  Technology, 
Durham,  New  Hampshire,  1991. 

[13]  J.K.  Hedrick,  V.  Garg,  J.C.  Gerdes,  D.B.  Maciuca,  and  D.  Swaroop, 

“Longitudinal  control  development  for  IVHS  fully  automated  and  semi-automated 
system  phase  III,”  Final  Report,  University  of  California,  Berkeley. 


53 


[14]  J.K.  Hedrick,  M.  Tomizuka,  and  P.  Varaiya,  “Control  issues  in  automated 
highway  systems,”  IEEE  Control  Systems  Magazine,  1994. 

[15]  D.V.A.H.G.  Swaroop,  “String  stability  of  interconnected  systems:  An  application 
to  platooning  in  automated  highway  systems,”  PhD  Dissertation,  University  of 
California,  Berkeley,  April  1994. 


54 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Infonnation  Center 
Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 


55 


